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ABSTRACT 


This thesis deals with the fault diagnosis of the 
hvdc system using pattern matching. The disturbances considered 
are ac voltage dips on rectifier as well as on inverter side and 
dc line faults occurring at different locations on the dc line. 
The system response to various disturbances has been obtained 
from the dynamic digital simulation of the hvdc system. The time 
domain discriminants have been calculated from these responses. 
The rectifier end direct current has been used to calculate these 
discriminants. Further, the log - transformed value of direct 
current has also been used to calculate the discriminants in 
order to improve the resolution in the values of the 
discriminants obtained for various disturbances. It has been 
found that the individual discriminants can not differentiate all 
of the disturbances. A combination of these discriminants have 
been used through a pattern matching technique to detect the 
various disturbances. The pattern matching technique, which has 
been used, is based on the concept of distance measure. 



CHAPTER 1 


INTRODUCTION 

The high voltage direct current (hvdc) 
transmission system has now emerged as a viable alternative to 
the high voltage alternating current transmission system for 
bulk power and long distance transmission and interconnection of 
ac systems. Uith the evergrowing use of the hvdc transmission 
system this has become the most promising area for the research 
oriented studies. This thesis is a step in this direction. 

One of the major advantages of the hvdc 
transmission is its fast controllability. The behaviour of an I 
hvdc system following a disturbance is predominantly affected by 
controller parameters. It has been proposed that the adaptive 
control for hvdc system shows better result than fixed parameter 
controllers [1]. However, the convergence time required for the 
adaptive control algorithm may be significant [2]. Further, it 
has been suggested that the look up table concept can overcome 
the problem of slow convergence associated with the adaptive 

I' 

control algorithm. One of the essential steps in the application | 
of look up table concept is fast state identification for making j 
decision about the choice of controller parameters. In the 
literatures [3-5], the system state identification using pattern I 
matching techniques have been suggested. It has also been 
reported that the pattern matching technique is computationally : 
time intensive. In this thesis an attempt has been made to i 
explore the possibilities of applying the pattern matching I 
technique for fault diagnosis of a hvdc system. The next section | 
deals with the basics of the pattern matching which has been S 




Fig. (lia): Learning before recognition 



FiQ- (1-lb): Learning and Recognition concurrently 
{Ret:''Computer oriented approches to pattern recognition“by W.S.Meisel) 















used subsequently for the fault diafinosis of a: hvdc system. 

1.1 PATTERN MATCHING 

Pattern matching begins with the class definition 
and labeled samples of those classes in some workable 
representation. The problem is solved when a decision rule is 
derived, which assigns a unique label to new pattern. There are 
two aspects of pattern matching namely developing a decision rule 
and using it [6]. 

The two stages of pattern matching, deriving a 
decision rule and using it , can be performed concurrently and 
sequentially. In Figure 1.1a the sequential procedure has been 
presented. In this method all the labeled pattern samples are 
collected and best decision rule based on those samples is 
derived. That decision rule is used without change to classify 
the unlabeled samples. Figure 1.1b shows the concurrent 
procedure. In this case the decision rule is modified as it is 
used. Here a sample pattern is presented and classified , an 
error detector indicates whether the classification. is correct or 
not and the decision rule is left unchanged or modified as 
appropriate. In this thesis the sequential procedure of pattern 
matching is used. 

The process, necessary in deriving the decision 
rule, is indicated diagramatically in Figure 1.2. The system from 
which the given pattern arises is characterised completely by its 
physical embodiment. Raw data describing the system is referred 
as the measurement space. The pattern space may be identical with 
the measurement space or several stages of intermediate 
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Fig. 1.2: Stages in the Derivation of decision Rule 




processing may be necessary. Feature selection (or preprocessing) 
is a process by which a sample pattern in the measurement space 
is described by a finite and usually smaller set of numbers 
called features. 

1.2 OBJECTIVE OF THE THESIS 

The need of a fast state identification technique 
has been established in the foregoing discussion. The main 
objective of the thesis is to develop a fault diagnosis technique 
using pattern matching to identify the various disturbances 
occurring in the two terminal high voltage direct current 
transmission system. 

1.3 ORGANISATION OF THE THESIS 

To meet the objective of fault diagnosis of the 
hvdc system the dynamic digital simulation of an hvdc system has 
been carried out for obtaining the system response to various 
disturbances. The digital computer representation of the hvdc 
system has been discussed in chapter 2. Both the conventional and 
unconventional control strategies of the hvdc system have been 
discussed. The recovery from the various disturbances has been 
validated with unconventional control at the inverter terminal. 

Chapter 3 deals with the time domain technique for 
the identification of disturbances. The various time domain 
discriminatnts have been discussed and their values have been 
obtained ' from the system response. The success of individual 
discriminants for the identification of disturbances has been 
examined. 

The pattern matching technique used for the 
identification of various disturbances has been discussed and 



utilized in chapter 4. Software for the fault diaenosia based on 
pattern ntatchine technique has been developed. The identification 
of the disturbances has been validated. Another system operatine 
condition (with unconventional control at inverter) has been 
considered to examine the effect of operatinfi condition on the 
identification process. 

The conclusions and future scope of this work have 
been presented in chapter 5. 



CHAPTER 2 


HVDC SYSTEM SIMULATION 


2.1 INTRODUCTION 

One of the major steps in the state identification 
procedure is to establish the knovledee about the system 
behaviour following a disturbance. Uith this objective the 
dynamic digital simulation of the hvdc system has been attempted. 
This chapter briefly describes the representation of the various 
subsystemsCthe ac system feeding the converters , converters and 
the associated controls and dc transmission network) for the 
purpose of simulation. In addition to the conventional control 
strategy of consteint current and extinction angle controls, the 
unconventional control (Constant Reactive Current control) for 
inverter has also been considered. Based on these control 
strategies the transient behaviour of the two terminal hvdc 
system has been studied to derive the appropriate knowledge for 
state identification. 

2.2 HVDC SYSTEM REPRESENTATION 

The approach employed for the representation of 
the hvdc system is to model each component or subsystem 
separately in a modular fashion. The various subsystem models are 
then interconnected using appropriate interfacing variables. The 
representation of various components of hvdc system is briefly 
described below. 

2.2.1 Converter Representation 

The converters which link two ac systems are 
represented as a variable voltage source behind a variable 



impedance. The converter equivalent circuit is shown in Fi£ure 
2.1. The voltage source (e^q) in equivalent circuit is a function 
of ac bus voltage and hence has to be calculated at each time 
instant. The equivalent circuit parameters R^q and L^q are 
dependent on the converter conduction pattern and are 
recalculated every time as the conduction pattern is changed. The 
source has been included to represent the effect of the dc 

system. Both the ac and dc voltage sources are the outputs of the 
ac and dc network models [ 12 ]. R^ and L 4 denote the resistance 
and inductance of the smoothing reactor. 

2.2.2 Control Representation 
(a) Conventional Control 

The hvdc converters are generally equipped with 
constant current and constant extinction angle controls. These 
are called as conventional controls. The controls of converters 
considered here are based on the digital technique proposed by 
Freris et al . [7]. The firing scheme is equidistant pulse control 
with pulse frequency control. The Inter Firing Period (IFF), 
which is the interval between two successive firing instant, is 
calculated as 

IFF = 60° + Qi 

where is the firing correction obtained from 

constant current or constant extinction angle control. Under 
steady state Qj^-0 and firing takes place at 60° interval. In case 
of constant current control, the firing correction is obtained as 

Ql » KVc 

where K is the gain and is the control signal. 
The block diagram for the calculation of IFF and firing 





correction is shown in Figure 2.2. The control signal nay he a 
linear function of current error as 

Vc = Kidd - Id ref) 

Arrillaga [8], however, has suggested that the 
inclusion of a derivative term in the calculation of control 
voltage improves the dynamic behaviour of the system. Then 
control signal becomes 

Vc = Kidd - Id ref) + K 2 (dld/dt) 

The control signal V^. is sampled and used for the 
calculation of firing correction. 

In case of constant extinction angle control the 
controller action takes place through two completely independent 
loops. These are inverter safety control and inverter optimum 
control. The extinction angle (F) is the current zero instant of 
a particular valve and the instant at which the voltage across it 
becomes positive again. The controller , by changing the firing 
instant , attempts to keep this angle at its specified minimum 
value. The change in the firing instant is accomplished according 
to the difference between the measured F and F reference. 

The inverter safety control acts when the measured 
F is less than the F reference and it can only reduce the firing 
angle. The firing correction is obtained as 

Qi = K3(F - Fref) 

The inverter optimum control tries to bring the 
system back to its optimum operating condition when the inverter 
safety is guaranteed. To obtain this, a record is kept of all the 
extinction angle measured in a cycle and if minimum of all these 
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INVERTER OPTIMUM CONTROL 


FIG. 2.3. CONSTANT EXTINCTION ANGLE CONTROL 





extinction angles is greater than T reference the firing angle is 
increased. Thus the inverter optimuin control operates only once 
in a cycle and the firing correction is obtained as 

Qi = 

The block diagram for the calculation of firing 
correction and IFF for constant extinction angle control is shovm 
in Figure 2.3. 

Cbl Unconventional Control 

To investigate the influence of other control 
scheme on the fault diagnosis , an unconventional control has 
also been considered for the inverter. The various unconventional 
controls are the power factor control, the constant reactive 
power control, the constant reactive current control etc. The 
unconventional control considered, is based on constant reactive 
current control suggested by Szechtman et al . [9]. 

Szechtman et al . [9] have shown that the constant 
reactive current control has a superior general behaviour 
compared to the conventional control. The inverter when in 
voltage control , which is the normal condition , exhibits a 
negative impedance characteristic that may affect the over all 
system stability. Weaker the ac system more pronounced is the 
negative impedance effect and therefore, the voltage instability 
problem is aggravated. The voltage instability problem could be 
minimized using compensation techniques. Thus the principle of 
minimizing reactive power consumption by operating inverter at 
minimum extinction angle is partially ineffective as it requires 
more voltage support. 

To circumvent the above problem , an alternate 



technique based on the cohtrol of reactive current at the 
inverter has been suggested «[ 9] . The firing angle of the inverter 
is regulated in a manner that the reactive current absorbed by 
the inverter remains ccxnstant. The calculation of firing 
correction and IFF is saine as in case of constant current 
control. The control signalc^^^ is generated as a function of 
converter reactive current: and reference reactive current. A 
derivative term has been' ^included to improve the dynamic 
behaviour of the system The control signal is obtained as 

Vc = KiCIq - + KjCdlq/dt) 

The values o'B/gain and K 2 have been chosen 
judiciously. The control signal is sampled and used to 
calculate the firing correction based on Pulse Frequency Control 
as 

Qi = KVc 

The firing c^Blrrection thus obtained is used to 
calculate the IFF which decides the next firing instant. Under 
steady state the firing corret:tion ($?£) is zero and firing takes 
place at dO^ interval. 

2.2.3 AC AND DC NETWORK REPRESENTATION 

The ac systemtls represented by an ideal voltage 
source behind a T - equivalent circuit as described in [10]. A 
three phase schematic repr^entation of the ac system and 
harmonic filters associated* with a particular terminal is shown 
in Figure 2.4. 

The inductance (L 3 ) and resistance (Rg) are 
determined at fundamental frequency from the knowledge of short 







circuit ratio and impedance angle at the converter terminal. The 
effective short circuit ratio gets modified in the presence of 
harmonic filters and/or shunt capacitors. The calculations of 
effective short circuit ratio and other ac system parameters have 
been reported in reference [11]. The shunt capacitor CCgD has 
been chosen based on reactive power requirement of the 
converters. The effect of converters on the ac system are 
represented by current sources and Ig. The tuned harmonic 
filters for 5 ,7 ,11 and, 13^ order harmonic with a second 
order high pass filter have been represented at the converter 
bus . 

The dc network comprises of dc filters , smoothing 
reactors and transmission lines. The transmission line is 
modelled as a Pi - equivalent circuit with the shunt arms 
consisting of capacitors, and series arm a combination of 
resistance and inductance. A more practical representation of 
transmission line is a number of Pi - section. This is done by 
dividing the transmission line into a number of sections, 
representing each section by its own Pi - network and connecting 
all such networks in series. In the dc system representation 
converters are represented as current source as shown in Figure 
2.5. The dc filters , used for filtering any harmonic present in 
the direct current, are also represented in dc system model. 

2 . 3 COMPUTER PROGRAM 

A computer program, based on the representation 
discussed in the preceding sections, has been developed and 
reported in reference [12]. This incorporates both 6 and 12 pulse 
converter representations along with analog firing control scheme 
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based on IPC and EPC. The program is extremely modular and 
therefore permits the implementation of different control schemes 
and other advanced features necessary to simulate the practical 
system. The program was augmented and reported in the reference 
[11], to incorporate the digital firing control scheme as given 
by Freris et al.[7]. The program is augmented to incorporate 
unconventional control (Constant Reactive Current Control) at the 
inverter terminal . The fault diagnosis subroutine has been 
developed, based on pattern matching technique, for the 
identification of various types of faults. The basic flow chart 
of augmented program is given in Figure 2.6. 

The simulation can start either assuming the 
system under steady state operating condition or with zero 
initial condition. In the former case the initial conditions for 
various state variables have to be supplied. The calculation of 
initial conditions for various state variables have been reported 
in reference [11]. 

2 . 4 CASE STUDY 

The augmented computer program is used to carry 
out various test simulations. The dynamic behaviour for various 
disturbances has been studied under unconventional control at the 
inverter terminal. The various types of disturbances given in 
Table 3,1 have been simulated. The recovery from these 
disturbances with constant reactive current control at inverter 
has been tested. The system response, for line to ground fault at 
inverter with conventional and unconventional, is shown in 
Figures 2.7 and 2.8. The recovery time with conventional and 
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unconventional control is same and it is approximately 280 ms. 
The recovery from three phase fault at inverter is shown in 
Figures 2.9 and 2.10 with conventional and unconventional 
controls respectively. 

2.5 DISCUSSION 

The digital computer representation of the hvdc 
system has been discussed in this chapter. Both the conventional 
and unconventional controls for the inverter terminal have been 
discussed. The digital computer program is augmented to 
incorporate an unconventional control ( Constant Reactive Current 
Control) at the inverter terminal. The recovery from the various 
disturbances has been studied. The system response to various 
disturbances, obtained from the dynamic digital simulation of the 
hvdc system with conventional and un conventional controls, has 
been used to obtain the time domain discriminants in the 


following chapters. 



CHAPTER 3 


TIME DOMAIN TECHNIQUE 

3.1 INTRODUCTION 

In this chapter various types of tine donain 
discr ininants have been discussed. These discrininants are used 
to identify various types of disturbances. The discrininants are 
calculated fron the system response to various disturbances 
obtained fron dynamic digital simulation of the hvdc system. The 
rectifier side direct current is used as a signal for the 
calculation of these discrininants. Another signal which is the 
log - transformed value of direct current has also been 
considered. The disturbances considered are the ac voltage dips 
and dc line faults occurring at different locations on the ac and 
dc sides. 

3.2 TIME DOMAIN DISCRIMINANTS 

The following tine donain discrininants have been 
considered for the identification of various disturbances [13]. 

(1) Mean 

The mean value of a signal is defined as 

Mean (p) =^--^ 

N 

where, = The value of signal at i instant 

N = The number of samples 
(2) Variance / Standard Deviation 

The second moment about the mean is called the 
variance of the signal. Mathematically the variance of a signal 
is defined as 



Variance (a) 


N 


5 . 


( Xi - M 


The square root of the variance is called the 
standard deviation of sifinal distribution. 

(3) Skewness 

Skewness, the third moment about the mean, gives a 

measure of skewness or asymmetry of the distribution. 

Hathematically skewness is defined as 

N 

^ M 

Skewness (S) = — -yr 

N v' 


(4) Kurtosis 

, Kurtosis, the fourth moment about the mean, gives 

the measure sharpness. Hathematically Kurtosis is defined as 

N 

1 iS, - “J* 

Kurtosis (K) = r 

N Q 


C5) Area/Integral 

^ Area under the distribution can be calculated by 

using Simpson’s formula for integration. The Simpson’s formula 

for integration can be approximated as 

' X h 

^ y(X) dX = — (Yq + 4Yi + Y2) 

J Xo 3 

Uhere , 

Xo = X(Xo) 

Yi = Y(Xo + h) 



Y2 = YCXo + 2h) 


h = Integration step length 
Y(X) = Distribution function 
This gives the area of the strip lying between Xq, 
X 2 and Yq , Y 2 - Area under the distribution can be obtained by 
summing the areas of such strips over the entire distribution. 
The total area has been normalized by dividing it with the number 
of samples being used. 

3.3 CHOICE OF SIGNAL 

The effect of a disturbance in a dc transmission 
network or in the associated ac system, is reflected on the dc 
line in general as a variation in direct current and direct 
voltage. Any of these two quantities can be considered for the 
calculation of various time domain discriminants. However, as the 
occurrence of the short circuit faults is more frequent than the 
open circuit faults. Thus the direct current forms a convenient 
basis for the calculation of various time domain discriminants. 
Following most of the short circuit faults the direct voltage 
will collapse to zero. This may cause considerable computational 
problem in the evaluation of the discriminants. Following a 
disturbance, the direct current is sampled and used for the 
calculation of the discriminants. The direct current samples can 
be obtained for a reasonable period of time, say, for about two 
cycles or so. 

3 . 4 CASE STUDY 

For the calculation of time domain discriminants , a 
detailed dynamic digital simulation of two terminal hvdc system 



Table 3 . 1 ; List of disturbances 


Case Discription 

No. 


1 . 
2 . 

3. 

4. 

5. 
6 . 
7 . 
8 . 
9. 

10 . 
11 . 
12 . 

13. 

14. 

15. 
16 . 


Single phase solid fault at inverter 
Single phase 80% dip at inverter 
Single phase 60% dip at inverter 
Single phase 40% dip at inverter 
Single phase 20% dip at inverter 
Two phase solid fault at inverter 
Three phase solid fault at inverter 
Remote three phase fault at inverter 
DC line fault in first Pi - section 
DC line fault in second Pi - section 
DC line fault in third Pi - section 
DC line fault in fourth Pi - section 
DC line fault in fifth Pi - section 
Single phase fault at rectifier 
Two phase fault at rectifier 
Three phase fault at rectifier 


has been carried out to obtain the tranaient reaponae of the 
direct current following various disturbances. The conventional 
controls for the converters have been considered here. The 
rectifier is equipped with constant current and ntinifflum alpha 
control. The inverter is equipped with constant current and 
constant extinction angle control. 

The various faults considered on the ac and dc 
sides are listed in Table 3.1. The direct current response for 
the various cases are given in Figures 3.1 to 3.16. The responses 
shown, are for a period of two ac cycles (40 ms) from the instant 
of initiation of fault. Utilizing the sampled values of the 
direct current, the time domain discriminants discussed in 
section 3.2, have been .evaluated and presented in Table 3.2. It 
can be observed from Table 3.2 that the values of discriminants 
are quite close for some of the disturbances. For example, the 
Mean value discriminant for two phase solid fault at inverter is 
1.375203 and that for single phase 80% dip at inverter it is 
1.374784. Similarly the value of Variance for single phase solid 
fault at inverter is 0.082278 and that for single phase 80% dip 
it is 0.082507. Thus due to the closeness in the values of 
discriminants for different faults, it may be difficult to 
differentiate the faults based on the individual discriminants 
alone. However, a combination of discriminants can successfully 
be used for identifying the type of fault. For an example, the 
two phase solid fault and single phase 80% dip, which may not be 
discriminated using Mean value, can be discriminated at the 
second level of decision based on the Variance. For some of the 
faults more number of decision levels may be required. This may 
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ultimately slow down the decision making process, and the 
advantage of using the time domain discriminants for 
identification may be lost. 

From the above discussion it is evident that in 
order to discriminate the various faults successfully either a 
combination of various discriminants should be used or if the 
identification is to be based on single discriminant alone then a 
choice of an alternate signal for the evaluation of discriminants 
should be explored. 

3.5 LOG TRANSFORMATION 

For calculating the various time domain 
discriminant, the log transformed value of direct current can be 
considered as an alternate signal. This can be obtained by 
performing the log transformation on the direct current. This has 
the following features [13]. 

(a) The direct current can be normalised to 1.0 
p.u. in the steady state and the logarithm of 1.0 is zero. Thus 
any definite value of log transformed direct current will 
indicate the disturbance in the system. 

(b) The data which is segregated can be compressed 
and wild variations can be smoothed out with the use of log 
transformation. 

(c) The log transformation performed on the signal 
of interest improves the resolution in the values of 
discriminants . 

3 . 6 CASE STUDY 

The log transformed value of direct current signal 
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Table 3.3: CContd. 
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has been used to calculate the time domain discriminants 
described in section 3.2. It can be observed from the Table 3.3 
that the sifin of the nean value indicates whether the fault is on 
rectifier side or on inverter side. The Hean value for two phase 
solid fault is 0.130203 and that for single phase 80t dip it is 
0.128440. Thus the Mean value obtained from the log - transformed 
value of direct current can discriminate these cases. It has been 
shown in section 3,4 that these cases were difficult to 
discriminate using Mean value obtained from the direct current 
signal. Further, the value of Variance for single phase solid 
fault is 0.008624 and that for single phase 80% dip it is 
0.008633. Due to the closeness in the values, the Variance can 
not discriminate these disturbances. 

Thus the discriminants obtained from the log 
transformed value of direct current also , can not identify all 
of the faults when the decision is based on individual 
discriminants alone. 

3.7 CONCLUSIONS 

The various time domain discriminants have been 
discussed and calculated. The signals used for calculation were 
the direct current and its log transformed value. It has been 
concluded that all of the faults can not be identified using 
individual discriminants alone. By using the discriminants 
obtained from both the signals simultaneously a good decision 
about the type of disturbance could be reached. However, a large 
number of comparisons may be required to identify a fault when 
the decision is based on individual discriminants alone. This may 
slow down the decision making process. 



CHAPTER 4 


PATTERN MATCHING FOR FAULT DIAGNOSIS 

4.1 INTRODUCTION 

The time domain discriminants have been calculated 
from the signal of interest and its log - transformed value in 
the last chapter. The problems associated with the identification 
of faults, when individual discriminants are used separately, 
have already been discussed in the previous chapter. To 
circumvent these problems various discriminants can be used in an 
organised manner through a pattern matching technique. The 
pattern matching techniques based on these discriminants are 
given in the literatures [5-6]. In this chapter a pattern matching 
technique has been discussed and used for the fault diagnosis. 

4.2 PATTERN MATCHING TECHNIQUE 

The pattern matching technique which has been used 
for the fault diagnosis is based on the concept of distance 
measure. This technique has been discussed below. 

4.2.1 Distance Measure 

The distance is the crucial concept in pattern 
matching. Closer a point to another point, more similar are the 
patterns represented by those points [6]. The patterns 
corresponding to various disturbances can be described by a set 
of time domain discriminants. These patterns form a domain. In 
this domain various disturbances are represented by points. The 
new pattern corresponding to an unknown disturbance, can also be 
represented as a point in the domain of patterns. The distances, 
between the point corresponding to new pattern and the points 



corresponding to the patterns available, are obtained. The new 
pattern will be similar to the pattern for which the distance 
obtained is minimum. 

The distance can be obtained using any distance 
function which satisfies the following constraints. 

(a) d(x,x) = 0 

(b) d(x,y) > 0 

(c) d(x,y) = d(y,x) and, 

(d) d(x,y) + d(y,z) > d(x,z) 

The conventional Euclidean distance function satisfies the above 
constraints. The Euclidean distance function can be given as 


dlk = 


where , 


r n 


Ilcxi - yik)^ 

1=1 



XI = The value of i^^ discriminant describing the 
new pattern 

’t#<h 

Vik “ The value of i^ discriminant describing the 
k pattern 

n = Number of discriminants 

1^11 

djk “ The distance between new pattern and the k 
pattern 

In the literatures [5-6], normalized samples of 
the patterns are used for obtaining the distance. It has been 
shown in reference [13] that the values of the time domain 
discriminants do not change significantly for samples/cycle more 
than 32. The two cycles information, from the initiation of 
fault, have been used to calculate these discriminants and hence 
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are used to describe a pattern then the total number of data to 
be stored will be 64. n, where, n is the total number of patterns. 
However, if the information contained in the samples of pattern 
is stored as a set of time domain discriminants then number of 
data to be stored will be considerably reduced. For an example , 
if three discriminants are used to describe a pattern the total 
number of data to be stored will be 3.n for n patterns. Hence, 
larfie number of memory space is required when the sampled value 
is used to describe the pattern. Further , this may slow down the 
decision makinfi process. Thus instead of describins a pattern by 
its sampled value it can be described by the discriminants 
calculated from these samples. 

4.3 THE SELECTION OF DISCRIMINANTS 

The pattern corresponding to various disturbances 
are described by a set of time domain discriminants. The elements 
in the set of time domain discriminants are Mean, Variance, 
Skewness, Kurtosis, Standard Deviation, Normalized Area and their 
log - transformed values. It has been shown in chapter 3 that the 
combination of these discriminants should be used to discriminate 
all the faults. If all of these discriminants are used 
independently, a large number of comparisons may be required in 
making decision about the type of fault. However, these 
discriminants when used through a pattern matching technique can 
successfully discriminate all the faults and the comparisons 
required may be less. The computational time, in making the 
decision, depends upon the number of discriminants in a set and 
the number of decision levels required to reach the decision. The 



number o£ the decision level depends upon how efficiently the 
discriminants have been selected. The number of discriminants in 
a set should be such that a right decision can be made with less 
number of decision levels. During the investigation it has been 
found that the minimum number of discriminants, which can 
successfully describe the pattern, are three. 

4.4 DECISION MAKING PROCESS 

To reduce the number of comparisons required to 
identify the type of disturbance, various Classes of disturbances 
have been defined. The Class may contain different Subclasses 
corresponding to various disturbances. It has been shown in 
chapter 3 that log - transformed value of direct current shows 
better trending. This property has been exploited in classifying 
the various disturbances. It has been found that the Mean, 
Variance and Standard Deviation obtained from the log 
transformed value of direct current could be used to describe the 
various Classes of disturbances. The various Subclasses are 
described by simple values of Skewness, Variance and Standard 
Deviation. The various Classes and Subclasses are given in Table 
4.1. This classification is based on the values of discriminants 
and has been developed through investigations. 

At the most, the two decision levels are required 
to identify a fault as observed from Table 4.1. For an 
example, if in the first level of decision it has been decided 
that the new pattern belongs to Class 6, the second level of 
decision will give the Subclass to which it belongs. Thus in two 
decision levels the type of disturbance can be identified. 



Table 4 . 1 : Clasaif ication of Diaturbancea 


Clasa 


1 

2 

3 

4 


5 

6 


7 

8 
9 


Subclaaa 


Description 


a 

b 


a 

b 

c 


a 

b 

c 

d 

e 


Single phase fault at inverter 

- Solid fault 

- 80% dip 

Two phase fault at inverter 

Three phase fault at inverter 

Remote single phase fault at inverter 

- 60% dip 

- 40% dip 

- 20% dip 

Remote three phase fault at inverter 
DC line fault 

- First Pi - section 

- Second Pi - section 

- Third Pi - section 

- Fourth Pi - section 

- Fifth Pi - section 

Single phase fault at rectifier 
Two phase fault at rectifier 
Three phase fault at rectifier 



4.5 RESULTS AND DISCUSSION 

A software has been developed which employs the 
pattern matchlnfi technique for the identification of faults. The 
Flow Chart for this has been given in Figure 4.1. The software 
developed has been incorporated in the existing digital computer 
program used for the dynamic digital simulation of the hvdc 
transmission system. The knowledge about the system behaviour 
corresponding to the various disturbances has been stored in 
terms of time domain discriminants. The various Classes of 
disturbances have been described by log - transformed Nean, 
Variance and Standard Deviation. In a particular Class, various 
Subclasses are described by simple Skewness, Variance and 
Standard Deviation. The reduced pattern space is given in Table 
4.2. 

The dynamic digital simulation of the hvdc 
transmission system including the software developed for fault 
diagnosis has been carried out. The identification of the various 

disturbances listed in Table 3.1 has been validated. The various j 

I 

test simulations has been carried out to examine the influence of i 
the fault duration, other control strategy and change in the ac 
system impedance on identification process. : 

I 

4.5.1 Influence of Fault Duration | 

Since, two cycles information about the direct i 
current variation has been used in the identification ; 
process, therefore, the fault diagnosis technique does not depend 
upon the duration of fault. The reduced pattern space given in 

j 

Table 4.2, has been used as knowledge about the disturbances. 
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READ FAULT No. (K) 

SET OF DISCRIMINANTS DESCRIBING IT 
total No of FAULTS CONSIDERED (N) 


calculate SET OF DISCRIMINANTS 
FOR THE UNKNOWN FAULT FROM THE 
SYSTEM RESPONSE 

INITIALIZE K = 1 


CALCULATE EUCLIDEAN DISTANCE 
d^K (x,y) 



diK(x,y) < MINDIST 


I FAULT = K 
MINDIST = 


K = K+1 


jYES 

WRITE I FAULT 


Flow Chart of Fault Diagnosis Software 











Table 4.2: Reduced Pattern Space (with conventional control 
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This has been ob'talned from the system response when the system 
is subjected to a disturbance of 3 cycle duration. The 
identification of the faults has been found to be independent of 
the duration of fault when the fault duration is increased to 4 
or 5 cycles- 

4.5.2 Influence of Control 

The proposed identification technique can be used 
to identify the disturbances when the system operating condition 
is changed. To illustrate this point, a different system 
operating condition has been chosen. The hvdc transmission 
systems are generally equipped with constant current and constant 
extinction angle controls. These are the conventional control 
strategies. The constant reactive current control, which is an 
unconventional control, has been considered for the inverter 
terminal. From the system response , with constant reactive 
current control for the inverter, the various discriminants have 
been calculated. The reduced pattern space has been given in 
Table 4.3. This has been used as a knowledge about the various 
disturbances. 

If the hvdc transmission system is equipped with 
conventional as well as unconventional control, the control 
strategy will always be known apriorl. However, the change in 
the control strategies may be detected as a fault. In either case 
the corresponding pattern space (data set) can be selected. The 
identification technique remains the same, only the knowledge has 
to be updated. 
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4.5.3 Influence of the Change in AC System Impedance 

The dynamic digital simulation of hvdc 
transmission system including the software developed for fault 


diagnosis , 

has 

been carried out with the 

change 

in ac 

system 

Impedance . 

It 

has been found that ilOt 

change 

in ac 

system 


impedance does not affect the system response greatly and various 
faults can be identified faithfully without change in the 
knowledge. If ac system impedance is changed such that it affects 
the system response to the various disturbances then also the 
various faults can be identified by selecting the appropriate 
pattern space. The fault diagnosis technique remains same. 

4.6 CONCLUSIONS 

A pattern matching technique has been discussed 
and used for the identification of various faults. The decision 
rule is based on the minimum distance concept. The dynamic 

[ 

digital simulation of the hvdc system including the software 

f 

1 

developed for fault diagnosis has been carried out. The following | 
conclusions have been made. | 

(1) The discriminants obtained from the log 
transformed value of signal can be used to define the various j 

I 

Claeses of patterns. I 

(2) Three discriminants are required to describe a 

pattern. 

i 

(3) The right decision about the type of fault can not 
be made until the knowledge, about the system behaviour stored as j 
a aet of discriminants, is correct. With the change in system | 

f 

operating conditions such as change in ac system impedance or i 



change in control strategy, the system response gets modified and 
corresponding pattern space should be selected. In any case the 
fault diagnosis technique remains same. 


CHAPTER 5 
CONCLUSIONS 

Thia thesis has mainly been devoted to the 
software development for fault diagnosis of the hvdc system based 
on pattern matching technique. 

The dynamic digital simulation of the hvdc system 
has been carried out to obtain the system response to various 
disturbances. The disturbances considered are the ac voltage dips 
and dc line faults occurring at different locations on the dc 
line. The system response thus obtained is used to calculate the 
time domain discriminants: The direct current pattern during the 
fault and its log - transformed value have been used to calculate 
these discriminants. The direct current response to various 
disturbances has been described by a set of discrimlnamts . A 
pattern matching technique based on distance measure has been 
used to identify the various faults. 

The hvdc transmission system is normally equipped 
with constant current and constant extinction angle control. 
These are the conventional controls for the hvdc transmission 
system. An unconventional control (Constant Reactive Current 
Control) for inverter termihAl has also been considered. This has 
been done to examine the influence of the system operating 
condition on the identification process. 

The salient findings of the present investigation 
ete given below: 

(1) The individual time domain discriminants can not 
differentiate all the faults. A good decision about the type of 






of discriminants 



individually. However, a large number of comparisons are required 
to identify the type of disturbance. This may slow down the 
decision making process. 

(2) The patterns of the direct current and its log 
transformed value corresponding to various disturbances can be 
described by a set of time domain discriminants, instead of its 
sampled value taken at appropriate Interval. This will reduce the 
memory requirement for storing the patterns. 

(3) A pattern matching technique based on the distance 
measure can successfully Identify the various disturbances. 

(4) The minimum number of discriminants which can 
describe ^ pattern successfully are three. 

(5) The fault diagnosis technique based on pattern 
matching is Independent of the duration of faults. 

(6) This technique can successfully identify the 
location of dc line fault. 

(7) Under the different system operating 
conditions , such as change in short circuit ratio of ac system and 
control strategy, for which the system behaviour following a 
disturbance gets modified a different set of knowledge 
corresponding to that operating condition is required. 

The major emphasis in the present study has been 
on the development of fault diagnosis technique using pattern 
matching. The various disturbances considered are the ac voltage 
dips on the rectifier as well as the inverter side and dc line 
faults. As an extension of this work in future, it would be worth 
studying the following aspects: 



(1) The identification of converter faults using pattern 

matching. 

(2) The decision rule can be made adaptive. 

(3) Development of an Expert System for the hvdc 
transmission system. 


T 
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APPENDIX A 


SYSTEM DATA (adapted from reference [ll] ) 

DC base voltage = lOO kV 
DC base current = 1 kA 

DC base impedance = 100 Ohms 

Rated Power of DC Link = 240 MW for 12 pulse system 

= 120 MW for 6 pulse system 

AC system frequency = 50 Hz 

Bridge Transformer: 

Resistance R^^ = R ^2 ~ Ohms 

Commutating Reactance = X ^2 ” 6.283 Ohms 

Smoothing Reactor: 

Resistance R^^ = ^^^12 ~ Ohms 

Reactance = X ^2 “ 314.159 Ohms 

Transmission Line: 

Total resistance = 8.64 Ohms 

Total inductance = 0.50148 H 

Total capacitance = 54.1645 pF 

Any number of PI sections upto a maximum of 10 may be spej 

i 

i, 

t 

Digital Controller - j 

a) Constant Current Controller I 

E 

= 0.4 radians/per unit current j 

^2 = 0.4 radians^/per unit current | 

b) Constant Extinction Angle Controller 

** 1 ® ■ ! 



0.25 


tlViS Voltages (L-M) 

= 42.78 kV, = 38.51 kV 

Valve turn off time = 3" 


a) 


b; 


Operating conditions: 
12 pulse: 


-d1 ' 'd2 

= 1000 

Amps 

a . = 5^ 

min 

"c = 

10‘> 


0 : 

f 

= 148.39® 

6 pulse; 

^d1 “ ^d2 

= 1000 

Amps 

a , =5% 

min 

"c = 

lOo 

= 12.88 

Q 

f 9 

= 148.38® 


DC Filter parameters; 
R = 0.096 p.u. 

= 0.20267 p.u. 

= 0.5947411 p.u. 

= 0.047183 p.u. 

K 0.136659 p.u. 

AC Filter Configuration; 




Series Tuned Filter 


High Pass Filter 



AC Filter Parameters 


Harmonic 

No. 


Resistance 

Ohms 


Inductance 
mH - 


5 1.11418 36.32860 

7 0.801103 18.16430 

11 0.504805 7,26573 

13 0.4225 5.20478 

HP 16.1428 3.00410 

VDCOL Parameters: 

6 Pulse Operation 
a) At Rectifier 


VDCOL Parameters: 
6 Pulse Operation 
a) At Rectifier 



= 0 . 00008 

sec J 

= 0.03 

sec 


Corner points: 





Voltage (p.u.): 

0.6 

0.2 

0.0 


Current (p.u.): 

1.0 

0.4 

0.4 

b} 

At Inverter 





* 0.00008 

sec f "^UP 

= 0.04 

sec 


Corner points: 





Voltage (p.u.): 

0.6 

0.2 

0.0 


Current (p.u.J: 

0.9 

0.3 

0.3 

12 

Pulse Operation 




a) 

At Rectifier 





* 0.00008 

sec; ^UP 

= 0.03 

sec 


Corner points: 





Voltage (p.u.): 

1 .2 

0.4 

0;0 


Current (p.u.): 

1 .0 

0.4 

0.4 



b) At Inverter 

" 0-00008 sec; T^,p = 0.04 

Corner points: 

Voltage (p.u.): 1.? 0.4 

Current (p.u.): 0,9 0,3 


sec 

0.0 

0.3 



1 0 1) u j 




date last 


D.teSI^0S905 

book is to be returned on the 
lamped. 
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